Observation of two-magnon bound states in the two-leg ladders of (Ca,La)i4Cu2404i 
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Phonon-assisted 2-magnon absorption is studied at T=4 K in the spin- 1/2 2- leg ladders of 
Cai4_iLai,Cu2404i (a::=5 and 4) for -E || c (legs) and -E || a (rungs). Two peaks at about 2140 
and 2800 cm~^ reflect van-Hove singularities in the density of states of the strongly dispersing 
2-niagnon singlet bound state, and a broad peak at « 4000 cm~^ is identifled with the 2-magnon 
continuum. Two different theoretical approaches (Jordan- Wigner fermions and perturbation theory) 
describe the data very well for Jy sslOSO — 1100 cm~^, Jy/ ~ 1 — 1.1. A striking similarity of the 
high-energy continuum absorption of the ladders and of the undoped high Tc cuprates is observed. 

PACS numbers: 74.72.Jt, 78.30.-j, 75.40.Gb, 75.10.Jm 



Low-dimensional quantum spin systems display a fas- 
cinating variety of low-energy excitations. Prominent 
examples are the fractional quantum states of one- 
dimensional (ID) chains, the S=l/2 spinous, or the va- 
riety of bound states in gapped spin liquids such as, e.g., 
dimerized chains or even-leg ladders. In 2D, the magnetic 
excitations of the undoped high- Tc cuprates are usually 
discussed in more conventional terms, namely renormal- 
ized spin waves. But these fail to describe the short- 
wavelength high-energy excitatioHS probed by optical 2- 
magnon-plus-phonon absorption.^ The nature of incoher- 
ent high-epergy excitations is currently under intensive 
debate.Erel 

In this light, the study of other cuprates with related 
topologies is very instructive. As far as high-energy ex- 
citations are concerned, the CU2O3 2-leg ladders realized 
in (Ca,La)i4Cu2404i (Ref. ||) provide a bridge between 
ID physics and the 2D Cu02 layers. Antiferromagnetic 
5 = 1/2 2-leg Heisenberg ladders are represented by the 
Hamiltonian 

W = { J|| (Si^iSi^i+i -I- 82,^82, ,;+i) + J_iSi^iS2,i} , (1) 



where J± and J|| denote the rung and leg couplings, re- 
spectively. For J||=0 one can excite local rung singlets to 
rung triplets which become dispersive on finite J|| . For 
J_L=0 the S—l chain excitations decay into asymptot- 
ically free 5=1/2 spinous. An intuitive picture of the 
"magnons" (elementary triplets) for J±, J|| ^ can be 
obtained from both limits: the elementary excitations 
are either dressed triplet excitations or pairs of bound 
spinous with a finite gap A as long as J_l > 0. In 
a gapped system it is particularly interesting whether 
bound states occur. Theoretical studies of 2-leg ladders 
show that bothpsinglet and triplet 2-magnon bound states 
always exist.ErtHl However, in the 5=1/2 copper oxides 



their experimental observation is a difficult task. Inelas- 
tic neutron scattering directly probes the spin gap, but 
cannot determine the high-energy excitations to a sut. 
ficient extent due to the large exchange interactions.tS 
Magnetic Raman scatteringEJ is restricted to k=0 exci- 
tations, but for the relevant values of J\\/J± the singlet 
bound state appears only for finite wave vector (see Fig. 3 
and Ref. [l^). Ift,this case, optical spectroscopy is the ap-, 
propriate tool.t2l Two-magnon-plus-phonon absorptiouEj 
is able to probe magnetic excitations throughout the en- 
tire Brillouin zone (BZ), since the simultaneous excita- 
tion of a phonon takes care of momentum conservation. 

In this paper, we present optical conductivity data 
of Cai4_a;Laa;Cu2404i and identify the 2-magnon singlet 
bound state by two peaks reflecting van-Hove singular- 
ities in the density of states (DOS) of the bound state. 
Our theoretical results for bound states in the experi- 
mentally relevant parameter range are based on two dif- 
ferent approaches, namely Jordan- Wigner fermionalj and 
perturbation expansion about the strong dimer couuliag 
limit up to 13th order using unitary transformations.tZlE£l 
Both yield an excellent description of the optical data for 
J||«1050 - 1100 cm-i and J||/J_l « 1 - 1.1. 

Single crystals with x—5 and 4 were—grown by the 
travelling solvent floating zone methodllB Their single 
phase structure and stoichiometry have been verified 
by x-ray, energy dispersive x-ray and thermogravimet- 
ric analyses.E3 The La content x determines the average 
oxidation state of Cu. A nominally undoped sample is 
obtained fO|M:=6, which probably is beyond the La solu- 
bility limit .t3 Single phase crystals could only be synthe- 
sized for X < 5. The samples with a;=5 and 4 on average 
contain n=l/24 and 2/24 holes per Cu, respectively. But 
x-ray absorption data show that at least for n < 4/24 the 
holes are located withinJJia. chains ,EJ which agrees with 
sEm Thus we consider the lad- 
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FIG. 1. Optical conductivity (T(tj) of Cai4_a;Lai:Cu24 04i 
(thick lines: a::=5; thin solid lines: x—4) for -E |j c (legs) 
and E \\ a (rungs) at 4 K. Dashed: Bimagnon-plus-phonon 
absorption of YBa2Cu306 is given for comparison/ Inset: 
(Ca,La)i4Cu24 04i data on a larger scale. Dotted: estimate 
of the electronic background (exponential fits for ui > 6500 
cm"^ for E\\a-uj> 4500 cm~^ for E\\c). 



ders to be undoped which is supported by the similarity 
of our results for x—5 and x=4. 

The optical conductivi1gi,(T(w) was determined by col- 
lecting both transmissionc^ and reflection data between 
500 and 12 000 cm~^ on a Fourier spectrometer. In Fig. |l| 
we display a{uj) of Cai4_a;Laa;Cu2404i at T=4 K for po- 
larization of the electrical field parallel to the legs, E \\ c, 
and to the rungs, E \\a. Phonon and multi-phonon ab- 
sorption dominates (t{lu) below « 1300 cm~^ (see inset). 
A steep increase of the electronic background, probably 
due to interband excitations of charge-transfer type, is 
observed above 6000 cm~^ for -E || a (3000 cm"-'^ for 
E II c). Our analysis focuses on the two peaks between 
2000 and 3000 cm'^ and the broad feature at 4000 cm"!. 
Also plotted in Fig. 1 is cr(w) of YBa2Cu306 (Ref. |l|), a 
typical example of the 2-magnon-plus-piionon absorption 
spectrum of the undoped 2D cupratesJlJ This compari- 
son gives a first motivation to interpret the peaks in (t(uj) 
of the ladders as magnetic excitations. Note that in the 
ladders both the exchange constants and the relevant Cu- 
O bond stretching phonon frequencies are comparable to 
those found in the 2D cuprates. Since the exchange cou- 
pling in the chains is much smaller, Jchain ~ — 14 cm^^ 
(Rcf . |2^) , their magnetic excitations do not contribute to 
<j{uj) in this frequency range. Reducing the La content x 
from 5 to 4 causes a slight blue-shift of the magnetic ab- 
sorption that is opposite to the red-shift of the electronic 
background. We attribute the blue-shift to an increase 
of the exchange constants caused by the reduction of the 
lattice parameters (Ca is smaller than La), and the red- 
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FIG. 2. Comparison of experiment (solid lines for x—5) 
and theory for J± — = 1050 cm~^ (OPT: optimized per- 
turbation) and 1100 cm~^ (JoWi: Jordan- Wigner fermions), 
with Uph ~ 600 cm^^. The exponential fits of the electronic 
background have been subtracted from the experimental data 
(dotted lines in Fig. |l|) . 



shift of the background to an increased hole density in 
the chains. 

Since spin is conserved, infrared absorption is sensitive 
to AS'=0, singlet excitations, e.g., the excitation of two 
5=1 magnons coupled to S'totmO. Direct excitation of 
two magnons is Raman activetj but infrared forbidden 
due to the inversion symmetry on a 180° Cu-O-Cu bond. 
We can effectively avoid this selection rule by simulta- 
neously exciting a CjitjO bond stretching phonon that 
breaks the symmetry.t£l Hence the lowest order infrared- 
active magnetic absorption is a 2-magnon-plus-phonon 
process. The low values of (j{uj) < 2 il^^cm^^ indeed 
indicate a weak higher-order absorption process. 

The 2-magnon-plus-phonon coutribution to cr(a;) was 
evaluated by Jurecka and BrenigEij in the extreme case 
of strong rung coupling J|| — s- 0. They found a single 
sharp peak reflecting the singlet bound state. In order 
to gain quantitative control over the experimentally rel- 
evant coupling regime of J\\/J± « 1 we use two different 
approaches to calculate the dynamical 4-point spin cor- 
relation function, naiaely Jordan- Wigner fermions and 
perturbation theory.EJ In the forpier, we make use of 
the Jordan- Wigner transformationE3 to rewrite the spins 
as fermions with a long-ranged phase factor. Expand- 
ing the phase factor yields new interaction terms be- 
tween the fermions which we treat in RPA (random 
phase approximation). The second approach is pertur- 
bative in nature. It|_isj)erformed by a continuous uni- 
tary transformation.EZrLj The 1-particle energies and the 
2-particle bound state energies are extrapolated by stan- 
dard techniques (Fade and Dlog-Pade) ; the spectral den- 
sities are computed by optimized perturbationEj Details 
will be given elsewhere. The magnon (elementary triplet) 
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FIG. 3. 1-Magnon dispersion (dotted), band edges of 
the 2-magnon continuum (solid), and bound singlet disper- 
sion (thick lines) for J||/Jx = l obtained from Jordan- Wigner 
fermions (left) and extrapolated perturbation (mid panel). 
Right: energy of the van-Hove singularities of the bound sin- 
glet at fc = TT and at the dispersion maximum as a function 
of J||/J±. 

dispersion, the 2-particle continuum and the dispersion 
of the singlet bound state obtained for J± = J|| are given 
in the two left panels of Fig. ^. The results of the two 
techniques agree well even though there is a quantita- 
tive discrepancy 10%) in some energies (dispersion 
at k=Q and tt, bound state energy). The magnon disper- 
sion of extrapolated perturhaiion (mid panel) reproduces 
previous results very welLcSEZl We focus on the singlet 
bound state which shows a considerable dip at fc = tt. 
This has gone unnoted thus far because it occurs only 
for J||/J_L ~ 0.5 (see right panel in Fig. 3 and Ref. pi] ). 

To compare with experiment we consider the effect of 
the phonons. The total spectral weight is obtained by 
taking into account a dependence of the exchange con- 
stants on the external electric field E and the displace- 
ments of O ions u, J||__L = J|| _|_(E, u) (Refs. p^,p]). The 
phonons modulate the intersite hopping and the on-site 
energies on both Cu and O sites. We expand J(E, u) to 
order (Pj/dudE which entails the coupling of a photon 
to a phonon and two neighboring spins. This determines 
how to integrate the spin response in the BZ. Here, the 
weight factor is a mixture of sin*(fc/2) (Ref. ^) and of an 
isotropic (Ljndependent) form factor, where the former 
dominates.t3 For simplicity we consider a pure form fac- 
tor sin*(A;/2) and Einstein phonons with ujph = 600 cm~^ 
as is common for the cuprates. Our findings depend only 
very weakly on the precise phonon energy. 

In Fig. H we compare the theoretical results with the 
experimental data. For the former an artificial broaden- 
ing of J±/20 is used. For the latter, the magnetic part 
of cr(uj) was obtained by subtracting the electronic back- 
ground (dotted lines in Fig. |l]). For E\\a, the background 
was determined unambiguously by an exponential fit for 
w > 6500 cm~^. The cra(^) curves for x—5 and 4 nearly 
coincide after subtraction of the background. This cor- 
roborates the assumption that the ladders are undoped. 
For E\\c, the measurable range of transmission was lim- 
ited to below « 4500 cm^^ by the higher absorption, and 



the precise shape and spectral weight of the highest peak 
at « 4000 cm~^ cannot be determined unambiguously. 

The two peaks at 2140 and 2800 cm^^ (-^==5) can be 
identified with the ID van-Hove singularities in the DOS 
of the singlet bound state. The lower peak corresponds 
to the singlet energy at the BZ boundary, E{k = tt), 
the higher one to the maximum of the singlet dispersion 
-Emax at about fc w 7r/2 (see Fig. ||). These two energies 
determine the two free magnetic parameters Jj_ and J|| . 
The experimental <t(w) is well described by both theories 
for J|, «1050 - 1100 cm-i and J||/Jj_«l (see Fig. 2) pro- 
viding an unambiguous identification of the experimental 
feature^ The quantitative analysis can be pushed a step 
furthered based on the values of i^max and E{k — tt) 
computed by extrapolated perturbation (right panel of 
Fig. 3). The strong dependence of i?max on J\\/J± allows 
to pinpoint the exchange constants, yielding the same 
value for J|| « 1050 cm~^ and a slightly larger ratio of 
J\\/ J± ~ l-l- Our interpretation of the experimental fea- 
tures is confirmed by the good agreement between theory 
and experiment concerning the line shape of the bound 
states in cr(uj). Excellent justification for this interpreta- 
tion is also provided by the selection rule stemming from 
reflection symmetry about the a axis (RSa). Both theo- 
ries show that the bound singlet at fc = tt is even under 
RSa. But the excitations at fc = tt are odd under RSa for 
E \\ a and even for -E || c. Thus the weight of the bound 
state varies as (fc — tt)^ for E \\a whereas it is prevailing 
for i? II c. This explains the low spectral weight of the 
lower peak for i? || a. It is reduced to a weak shoulder. 

Shortcomings of the theory are the overestimation of 
the spectral weight of the 2800 cm~^ peak for || a 
and that the onset of (j{u!) around 2000 cm~^ is sharper 
than observed experimentally. However, the agreement 
is better than one may have expected since we neglected 
both the frustrating coupling between neighboring lad- 
ders and the ring exchange. A finite inter-ladder cou- 
pling will produce a dispersion of the bound state along 
ka and thereby broaden the features in (j{w), which can 
explain a smearing out of the onset at 2000 cm~^. With- 
out the ring exchange, the analysis of experimental data 
of various other techniques suggested a larger leg cou- 
pling J\\/J± ^ 1.5 (for a detailed discussion, see Ref. 
p^ . One finds a 1-magnon gap of A w 280 cm^^ and a 
dispersion extending up to w 1550 cm~^ (Ref. |l^). For 
J\\/J±—l the dispersion extends from A « 0.5J_l up to 
~ 2J|| (see Fig. ||). One obvious way to reduce A with 
respect to the maximum is to increase the ratio J\\/Jj_, 
in this case to ^ 1.5. Such a large value is difficult to 
reconcile with the microscopic parameters, in particular 
with the similar Cu-0 bond lengths along the legs aud. 
the rung which has provoked a controversial discussion.cj 
Our analysis shows that any values of larger than 

1.2 can be excluded. Recently, it was pointed out that 
the neutron data are also consistent with an isotropic ex- 
change J\\/J± ~ 1 — 1.1 and J|| w 900 cm~L.if a ring 
exchange of « 0.15 J|| is taken into account Jl3 Consid- 
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ering the fact that the ring exchange will renormalize J|| 
and Jl, this is in perfect agreement with our findings. 

Judging the weight of the high-energy continuum one 
should note that it increases by taking into account an 
admixture of an isotropic form factor to the sin"'(fc/2) 
form factor used here. Second, the 2-particle excitations 
considered here carry only « 75% of the totaLspectral 
weight. In the 2D cuprates, Lorenzana et aLtiJ argued 
that ring exchange increases the spectral weight at high 
energies. It will be most interesting to check the influ- 
ence of ring exchange on the continuum of the ladders, 
since these provide a technically much better controlled 
ground. Note that the theoretical curves in Fig. |2| de- 
scribe the frequency range of continuum absorption very 
well, which gives strong support to our interpretation. 

Finally, we address the similarity oi cr{uj) of the ladder 
(Ca,La)i4Cu2404i and of the undoped 2D cuprates (see 
Fig. 0). The 2D case does not show a truly bound state, 
the sharp peak at 2800 cm^^ in a{uj) of YBa2Cu306 
is caused by a resonaiice, an almost bound state lying 
within the continuumJij This main peak is well described 
in terms of 2-magnon-plus-phonon absorption, but the 
high-euergy peak at 3800 cm~^ is absent in spin- wave 
theory.El Note that this discrepancy is particular for the 
2D 5=1/2 case, the high-energy excitations are absent in 
the comparable S=l system La2Ni04 (Refs. [l5| , ^ ). In 
the cuprates, the magnetic origin of both peaks has been 
confirmed by abspsption and Raman measurements un- 
der high pressure.E3 It has been suggestedEJ that the high- 
energy weight in 2D 5=1/2 compounds is due to strong 
quantum fluctuations that go beyond spin-wave theory. 
The intriguing similarity of the 2-particle continuum of 
the 5=1/2 quasi-lD ladder (Ca,La)i4Cu2404i with the 
3800 cm~^ peak of YBa2Cu306 strongly indicates that 
the high-energy spectral weight is indeed a signature of 
strong quantum fluctuations in this compound. 

In conclusion, the existence of the singlet bound state 
and the 2-particle continuum was demonstrated exper- 
imentally in the optical conductivity spectrum cr(w) of 
the 5=1/2 2-leg ladder compound (Ca,La)i4Cu2404i. 
By two independent theoretical approaches we confirmed 
that the two sharp peaks plus the important broad con- 
tinuum seen in experiment are the unambiguous signa- 
ture of a ID dispersive bound singlet with strong inco- 
herent quantum fluctuations. The experimental spectral 
weight distribution reflects the theoretical selection rules 
perfectly. Quantitative analysis yields Jy « 1050— 1100 
cm~^ and J\\/J± « 1 — 1.1. Our findings indicate that 
the similar experimental results for the high-energy exci- 
tations of the undoped high- Tc materials are also due to 
strong quantum fluctuations. 
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